Unusual Zeeman Field Effects in 2D Spin-Orbit Coupled Fermi Superfluids 
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We investigate the Zeeman field effects on the bulk superfluid properties and the collective modes in 2D at- 
tractive atomic Fermi gases with Rashba-type spin-orbit coupling. In the presence of a large spin-orbit coupling, 
the system undergoes a quantum phase transition to a topological superfluid state at a critical Zeeman field. We 
show that the non-analyticities of the thermodynamic functions as well as other physical quantities at the quan- 
tum phase transition originate from the infrared singularities caused by the gapless fermionic spectrum. The 
same argument applies also to the BCS-BEC evolution in 2D fermionic superfluids with p- or d-wave pairing. 
The superfluid density «j and the velocity of the Goldstone sound mode behave oppositely in the normal 
and the topological superfluid phases: they are suppressed by the Zeeman field in the normal superfluid phase, 
but get enhanced in the topological superfluid phase. The velocity of the Goldstone sound mode also shows 
non-analyticity at the quantum phase transition. For large Zeeman field, we find ;!j — > n and c, — > vp, where n 
is the total fermion density and up is the Fermi velocity of non-interacting system. The unusual behaviors of the 
superfluid density and the collective modes can be understood by the fact that the spin-orbit coupled superfluid 
state at large Zeeman field can be mapped to the p, + ipy superfluid state of spinless fermions. 
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The Zeeman field (ZF, denoted by h) effects on the BCS 
superconductivity have been a longstanding problem for sev- 
eral decades At weak coupling, the BCS state undergoes 
a first-order phase transition to the normal state at hcc - 
0.707Ao where Aq is the pairing gap at h - 0. Fur- 
ther studies showed that the inhomogeneous Fulde-Ferrell- 
Larkin-Ovchinnikov state l3j| survives in a narrow window 
between hcc and /zfflo = 0.754Ao. The ZF effects on the 
fermionic superfluidity in the whole BCS-BEC crossover jUj] 
regime have been experimentally studied in recent years jsl]. 
Two-component atomic Fermi gases with population imbal- 
ance were realized to simulate the ZF effects. Around the 
Feshbach resonance, the phase separation between the super- 
fluid and the normal phases has been observed in accordance 
with the first-order phase transition. Despite the rich phase 
structure in the BCS-BEC crossover the superfluidity is 
completely destroyed at large enough ZF. 

Recent progresses on synthetic spin-orbit coupling (SOC) 
for neutral atoms iQ-lst] provide new way to study the SOC 
effects on the fermionic superfluidity llioll . Previous stud- 
ies of 2D solid-state systems showed that the SOC induces 
spin-triplet pairing, even though the attractive interaction is s- 
wave 111 ill . By applying a large ZF, the 2D system undergoes a 
topological phase transition to a topological superconducting 
state, where the non-Abelian topological order and Majorana 
fermionic modes can be realized 11211 . However, the proper- 
ties of the bulk phase transition and the collective modes are 
less understood for such systems. 

In this paper, we study the bulk phase transition and the 
collective modes in 2D atomic Fermi gases with combined 
SOC and ZF effects. The main results can be summarized as 
follows: (i) The bulk phase transition originates from the in- 
frared singularities caused by the gapless fermionic spectrum. 
The analyticity of any physical quantity across the phase tran- 



sition can be determined by analyzing the infrared behavior of 
the momentum integrals. For the present system, we find that 
the quantum phase transition is of third order (ii) The super- 
fluid density « , and the velocity of the Goldstone sound mode 
Cs behave oppositely in the normal and the topological super- 
fluid phases. They are suppressed by the ZF in the normal 
superfluid phase, but turn to increase with the ZF in the topo- 
logical superfluid phase. The sound velocity c, also shows 
non-analyticity across the phase transition, (iii) For very large 
ZF (h oo), we obtain analytically «j — > n and Cs — > vp, 
where n is the total fermion density and vp is the Fermi ve- 
locity of non-interacting systems. We show that the unusual 
behaviors of the superfluid density and the collective modes 
manifest the fact that the spin-orbit coupled superfluid state 
at large ZF can be mapped to the + ipy superfluid state of 
spinless fermions. 

Model and effective potential - The many-body Hamilto- 
nian for the 2D Fermi system we considered can be written as 
H = Hs + //int, where 

//, = J d~rif,\r) J|- - ^ + + 'HzJ ifr(r), 

//in, j ^'\{v),p\{v)if,^{v)if,^(r). (1) 

Here i/'(r) = [i^-f(r), i/rj^(r)]^ represents the two-component 
fermion fields, p = p ^e , + ^,6,, is the 2D momentum oper- 
ator with Pi - —ittdi, (T - cTfe V + cr,,ev with cr,- being the Pauli 
matrices, and ji is the chemical potential. The contact cou- 
pling U > Q denotes the attractive s-wave interaction between 
unlike spins. The ZF term reads 'Hz - -hcr^ and the spin- 
dependent term 'Hso - Act ■ p is the 2D SOC 113]. We set 
h > and A > without loss of generality. In the following 
we use the units Ti - - M - I. 
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In the imaginary-time functional integral formalism (tem- 
perature T - IIP), the partition function of the sys- 
tem is J3 = exp |->S[i^, (/'']} with the action 

S[il/,if/''^] = jf^fT[/^f2ri/r"l'5ri/' + //( (/',</'"'")]■ Introducing the 
pair field <b{x) - -Ui//i{x)tlr-i{x) [x - (r, r)] and integrat- 
ing out the fermionic degrees of freedom, we obtain X = 
JSOSO' exp{ - >Seif[<[>,<[)' ]}, where the effective action is 
given by 

S,d^, O' ] = -i J dxmx)\- - lTrln[G-'(x, x')]. (2) 

In the Nambu-Gor'kov representation, the inverse single- 
particle Green's function reads 



G-\x,x') = 



G-\x) 
<[)"l"(x) 



0{x) 
G-J(x) 



6(x — x'). 



(3) 



where G^^{x) - -dj + licr^ + (p-/2 + Acr ■ p - /d). 

In the superfluid state, the pairing field <l)(x) acquires a 
nonzero expectation value (<l>(x)) = A which we set to be 
real without loss of generality. By separating the pairing 
field as <l>{x) - A + 4>{x), the effective action >Seff[0, 0' ] 
can be expanded in powers of the complex fluctuation field 
We have >Seff[<l), <D^] = >Sfjf'(A) + >Sfj^[0, 0"'"] + ■ ■ ■ , where 
>S*^!(A) = >Seff[A, A] is the saddle-point or mean-field effective 
action with the pair potential A determined by the saddle point 
condition dS^^^jdA - 0. The collective modes are determined 
by the Gaussian-fluctuation part S^^^icj), 0^]. 

Infrared singularity and bulk phase transition - The single- 
particle excitation spectra can be read from the pole of the 
fermion Green's function 0(K), which is obtained from G by 
the replacement O ^ A. Here K = {ico,,, k) with oj,, being the 
fermion Matsubara frequency. Working out the explicit form 
of Q{K), we obtain the quasiparticle dispersion {a - +), 
where is^ is given by E"^ - [E^ + '7k + 2Q'^k]'^^- Here we 
have defined £k = [^l + A^]^'^, Tjk = [^^k^ + h'^^l^, and 
^k = [^k'?k + h^^^]^'^ with = kV2 - n. From the identity 
{E^f(E:^f = {El-T]lf+AA^\i^A^, we find that the fermionic 
excitations are fully gapped for A + Q except the case that the 
condition Co = /U^ + A- - = is satisfied. For Co = 0, the 
lower branch has a linear dispersion near k = 0, that is, 
E:^ - Wclkl + (9(|kp), where the velocity Uc - AA/h. 

The gapless fermionic spectrum causes non-analyticities of 
some physical quantities at the critical point Co = 0. To be 
specific, we consider the thermodynamic potential Q.(fi, h) = 
Q.{fi, h, A(/i, h)) at zero temperature, where 



^k+^k 



k^ + fB 



+ ^k 



(4) 



Here we have used the usual regularization ?7 ' = I]k(k^ + 
ffi) ' for 2D systems with being the binding energy of 
the two-body bound state in the absence of SOC. To obtain the 
thermodynamic potential h), the pair potential A(yU, h), 
which is regarded as an implicit function of fi and h, should 
be determined by the gap equation dQ.{ii, h, A)/5A - 0. 



To study the analyticity of the thermodynamic poten- 
tial or its derivatives with respect to ji and h, we con- 
sider two susceptibilities - -d^Q.(jj.,h)/dfi^ and Xhh - 
are related to the isothermal compress- 
ibility and the spin susceptibility, respectively. To obtain their 
explicit expressions, we need the derivatives dA(p, h)/dfj. and 
dA(p, h)/dh. They can be obtained from the gap equation 
dQ.{jj,h,A)/dA = 0. Finally, the two susceptibilities can be 
evaluated as 



Xhh 



dn(p., h. A) 1 / dn(p., h, A) 

dii dA 
dmQi,h,A) 1 (dm(fx,h,A) 
A 



dh 



dA 



(5) 



Here A = d^D.(^i, h, A)/dA^, n = -5Q(jU, /i, A)/5yU is the total 
density, and m = -dQ.{fi, It, A)/dh is the spin polarization. 

We find that the expressions of Xf^it and Xhh contain some 
momentum integrals of the following type 



Jo 



kdk-^g{k). 



(6) 



where (3i = 1 - /i-/^k, Qi ^ I - /?^/^k, and ^ I - eI/^u- 
The function g(k) approaches some nonzero constant for k — > 
0. At Co = 0, the integrals J,y are infrared safe since the 
quantities (3, go as k^ for k ^ 0. Therefore, Xt^ti and Xhh 
are continuous across the phase transition. However, the Z-th 
derivatives of the susceptibilities with respect to fi or li contain 
momentum integrals whose infrared behaviors go as 



ni; jA-2l p 



dkk 



2-21 



(7) 



For I = 2, the infrared divergence shows up. Therefore, the 
fourth derivative of Q(//, li) is divergent at the phase transition. 
Then the third derivative is discontinuous and hence the sus- 
ceptibilities show non-analyticities. Based on these observa- 
tions, we conclude that the quantum phase transition at Co = 
is of third order 11511 . 

For homogeneous system, the pair potential A and the 
chemical potential /i are determined by imposing the total 
density n - k^l{2n) — ep/n. The system can be character- 
ized by two dimensionless parameters: the attractive strength 
ln(^Ffl2D) and the SOC strength A/kp. Here the 2D scatter- 
ing length a2D is defined as - '^e^^'' li.Ma^^ fl^ with 
y - 0.577216 being the Euler's constant. The numerical 
results presented in this paper are for ln(A:Ffl2D) = 2 and 
Ajkp - 0.5. The quantum phase transition occurs at h - he - 
-^yU^ + A- ^ 0.5 lep. Increasing the attraction and/or SOC en- 
hances the pairing potential and hence the critical field he, but 
does not lead to qualitatively different results. As shown in 
Fig. [T] the pair potential A, while is suppressed by the ZF, goes 
smoothly but never vanishes at large h. The chemical potential 
H goes smoothly and reaches a maximum at the phase transi- 
tion. For h > he, the system is a topological superfluid 11211 . 
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FIG. 1: The pair potential A (a), the chemical potential \i (b), the sus- 
ceptibilities^^^ and^;,;, (c) and the bulk excitation gap Eg (d) as func- 
tions of h. The dashed lines denote the critical ZF = -sjff + A^. 



Fig. [Tfc) shows the susceptibilities andxui- They are con- 
tinuous but not smooth at the phase transition, as we expected. 
Fig- Illd) shows the bulk excitation gap Eg = mmj^lE^, E^}. 
It equals the pair potential A only at h - 0. Near the phase 
transition, it goes non-monotonously. 

Superfluid density and collective modes - To study the 
behaviors of the superfluid density and the collective 
modes across the quantum phase transition, we consider the 
Gaussian-fluctuation part >S*^![0, 0^]. It can be written in 
a bilinear form 5 A'(e)M(0A(g), where Q = {iv„,q) 
with v„ being the boson Matsubara frequency, h{Q) - 
[0(G), 4>^i-Q)]^, and the 2 X 2 matrix M(0 is the inverse of 
the collective-mode propagator. The matrix elements of M(2) 
are constructed by using the fermion propagator 0(K). We 
haveM„(0 = M22(-e) = ^ + i Tr [^„(/: + 0fe(/:)], 
Mi2(0 = ^Y^K^^mxiiK + Q)QniK)l and M2i(0 = 
5 IjK Tr [^21 + 0^21 (^)]- Taking the analytical continu- 
ation iv„ — > w H- iO^, the dispersions w(q) of the collective 
modes are determined by the equation det M[(y(q), q] = 0. 

We can decompose Mn(w,q)asMn(w,q) = M|j(dLt,q)-)- 
Mj"j(w, q), where M||(w, q) and Mj^|(w, q) are even and 
odd functions of <jj, respectively. Meanwhile Mi2(<y, q) and 
M2i(<w, q) are even functions of oj and can be expressed as 
Mi2((y,q) = M*j(w,q) = M\^(u),q) + iM^^(oj,q). The term 
Mj2((y, q) oc hA^ vanishes when h or A is zero. Then we de- 
compose the complex field (p{x) into its amplitude mode p{x) 
and phase mode 0{x), <p(x) - p(x) + iAd(x). The effective 

I '2') 

action S J then takes the form 

err 



iZ(p(-e)^(-e))N(0(2jg 



(8) 



2(M| 



+ Ml.), 



where the matrix N(0 reads Nii(0 - ^vir^jj 1 ^,^^2 
N22(0 = 2A2(M+ - MI2), Ni2(e) = 2/A(M-; - M-^), 
and N2i(g) = -2/A(Mj"j + iM-^). Since Mj"j(0,q) = and 
Mj^2(dLi, 0) - 0, the amplitude and phase modes decouple com- 
pletely at (o), q) - (0, 0). At the saddle point we have precisely 
M+ (0,0) = M|2(0,0). Therefore the phase mode at q = is 
gapless, that is, the Goldstone sound mode or the Anderson- 
Bogoliubov mode for neutral Fermi superfluids. 

To study the low energy behaviors of the collective modes, 
we make a small q and o) expansion of N((2) at zero tempera- 
ture. In general, the expansion takes the form Ni 1 - A + Cq^ - 
Dci? + - ■ ■ , N22 - Jq--RcL)^-\ — , andNi2 - Njj - -iBa)+- ■ ■ . 
The term Mp(w, q) does not contribute in this expansion. The 
explicit forms of the expansion parameters are given by 

\2 



B = - 



V V 



k 

4^k 



{Elf 



\ +a- 



1 + ff- 



A\^ 



^k^k- 



El 



^k 



pa 

^k 



D 



■zz 

a=± k 



(Elf - A- 



A^k^^. 



A^ A^k}h^ 



{Elf 



k 



R 



k ^K^Kf a 



(I 



1 + 



^■k 



^k^k 



y t 

V(^k+^k)' 



1 + 



El - tjI 2A\^A^ 



^k 

A- A\^^', 



KE-k 



n 
4M 

+ a 



(Elf d 



—yy— 



1 + 



h-El A^k^h^A^ 



£■2 1 

^k 



2^k 



(9) 



The parameter A equals the quantity d'H/dA- in (5) at the 
saddle point. The phase stiffness J is related to the super- 
fluid density by 7 = n J{4M) {M = 1 in our units). «s 
can also be obtained from its standard definition UtIi . When 
the superfluid moves with a uniform velocity y,, the pair field 
transforms as ^ (j)g2/MV, r -pj^g superfluid density n , is de- 
fined as the response of the thermodynamic potential f2 to an 
infinitesimal velocity u^, that is, ^(i/j) = n{0)+\nsV^+O{v'l). 

Analyzing the infrared behaviors of the momentum inte- 
grals, the analyticities of the expansion parameters across the 
phase transition can be summarized as follows: (1) The phase 
stiffness J and hence the superfluid density is smooth; (2) 
The parameters A,B,R are continuous but not smooth; (3) D 
is divergent. The numerical results for these expansion pa- 
rameters and the sound velocity Cj = -\Jj/(R + B-/A) in the 
homogeneous system are shown in Figs. [2] and [3] Note that 
the superfluid density does not equal the total density n even 
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FIG. 2; The expansion parameters A, B, D and R as functions of h. 



at /i =0 due to the lack of Galilean invariance in the pres- 
ence ofsoc iH. Due to the non-analyticities of A, B and R, 
the sound velocity also shows non-analyticity at the phase 
transition. Moreover, we find that n.s and Cs behave oppositely 
in the normal and the topological superfluid phases. They are 
suppressed by the ZF in the normal superfluid phase, but get 
enhanced by the ZF in the topological superfluid phase. This 
is quite unusual since we generally expect that the superflu- 
idity should be suppressed by the ZF. On the other hand, the 
divergence of D indicates that the amplitude or Higgs mode 
becomes a soft mode around the phase transition. 

Analytical results for large Zeeman field - To understand 
the unusual behaviors of and Cs in the topological super- 
fluid phase, it is useful to reexpress the mean-field theory in 
the helicity representation ifisll . The helicity basis 
is related to the ordinary basis (i/ff, i/rj)^ by a SU(2) trans- 
formation. In the helicity basis the single-particle Hamilto- 
nian is diagonal, that is, Hs - 2„=± 2k^k'^ff(k)iA(!'(k) where 
^jj^ = + Q''7k- Therefore, the system can be viewed as a two- 
band system. The ZF provides a band gap 2h at k = 0. In 
the presence of pairing, the mean-field approximation for //int 
reads //int - \ T,a,i3=± 2k[AQ./;(k)i/r;Uk)iA^(-k)+/Lc.]. The new 
k-dependent pair potentials Aapi^i.) read A+_(k) = -A_+(k) = 
-As(k) and A++(k) - A*_(k) - -At(k), where the interband 
and the intraband pair potentials are given by As(k) - hA/r]]i 
and At(k) - Aik^ - iky)A/T]k. Using these new pair potentials, 
the quasiparticle dispersions can be expressed as 

£k = 7([^t + |A,(k)P]l/2 + ;,k)2 + |A,(k)p. (10) 

For h » he, we find that the pair potential goes as A ^ 
a/li^, while the chemical potential ^ -h + b, where a and 
b are some constants and b «: Therefore, the upper band 
with dispersion has a large gap and essentially plays no 



FIG. 3: The superfluid density 77, (in units of n) and the velocity of 
the Goldstone mode c, (in units of the Fermi velocity up = kp/M) as 
functions of h. 



role in fermion pairing. The lower band opens a Fermi 

surface atk^kp = [2(A^ -(-//+ ^JA* + 2A^-^ + li^)]^'^. Since 
the pair potential A <K h, the total density n is carried by the 
lower band. We have n - 2k ®i^k - ^k) - k^/i'^TT) and hence 
- V2^F- Then the system can be regarded as a weakly 
coupled + ipy superfluid of spinless fermions where the 
pairing occurs around the Fermi surface k = kp. The interband 
pair potential As(k) can be safely dropped and we have ^ 
[(^k)^ + |At(k)p]'/^. Near the Fermi surface, we get E^ ^ 
[Vp{k - kp)^ + E^] '^^ where the Fermi velocity vf - V2l'f( 1 - 
A'/rip) and the bulk excitation gap reads Eg ^ AAkp/rjp. Here 
we have defined rjp - [A^k^ + /i^]'^^. 

Based on the above observations, the superfluid density Hj 
can be approximated as - n(l - A~ jrip). Therefore, for 
h — > oo, we have n., — > «. It manifests the fact that, for large h, 
the pairing occurs only in the lower band which carries nearly 
the total density. Meanwhile, the other expansion parameters 
A, B and R are dominated by the terms that are peaked at the 
Fermi surface k - kp. Using the same integral technique in 
BCS theory, we obtain B^jA ^ and 7? ^ l/[8;r(l - A^jrip)]. 
Therefore, the sound velocity — » vp for h oo. This result 
can be reexpressed as Cj ^ vpj y/2, which is just the sound 
velocity of weakly coupled 2D Fermi superfluids. 

In summary, as the ZF is increased, the system behaves 
more and more like a p ^ + ipy superfluid of spinless fermions. 
Therefore, the fermion pairing in the topological superfluid 
phase feels less stress than in the normal superfluid phase. 
This explains the unusual behaviors of n , and c, at large ZF. 

Indication for p- and d-wave pairings - Finally, we 
point out that the infrared singularities which cause the non- 
analyticities should also show up in other systems, such as 
the 2D BCS-BEC evolution with p- and (i-wave pairings ll20ll . 
In such systems, the single-particle excitation spectrum is 
Ek = [^l + |A(k)|2]'/2, where A(k) ~ k for p-wave and 
A(k) ~ k^ for c/-wave pairings. At the quantum critical point 
/i = 0, the dispersion at low k goes as E\i ~ k for p-wave 
and Ek ~ k~ for tZ-wave pairings. Therefore, we expect that 
the collective-mode properties in such systems also show non- 
analyticities. The non-analytical behaviors of the collective 
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modes can be measured by using the Bragg spectroscopy 12111 . 
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